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RESUMEN 

Una muestra de 229618 galaxias con lmeas angostas de emision es utilizada 
para establecer dos nuevos tipos de evidencias inambiguas de agujeros ne- 
gros supcrmasivos en el centro de sus niicleos: 1) las galaxias Seyfcrt 2 y los 
LINERs siguen la misma Icy de potencia caracteristica relacionando la lumi- 
nosidad del flujo ionizante con la del continuo; 2) ambas muestran la mAjs alta 
concentration de masa en sus centro, independientemente de la morfologia de 
la galaxia, consistente con altas energias de ligation. El ancho a potencia 
media esta relacionado con la concentration dc masa, sugiriendo que la en- 
ergfa cinetica del gas en los AGNs tienc un origen gravitational. Dcntro del 
i-Ch ■ modelo estandar de acrecion, los Objetos de Transition, las galaxias Seyfert 

2 y los LINERs representan AGNs formando agujeros negros supcrmasivos a 
diferentes escalas en masa, o prodrian estar relacionados dentro de un proceso 
de evolution, donde los LINERs serfan el producto final. 

ABSTRACT 

A sample of 229618 narrow emission-line galaxies is used to establish two 
new unambiguous type of evidence for supcrmassivc black holes at the center 
of their nuclei: 1) the Seyfert 2 galaxies and LINERs follow the same char- 
acteristic power law relating the luminosity of ionized flux with that of the 
continuum; 2) both show the highest concentration of mass at their center, 
independent of the morphology of the galaxy, consistent with higher binding 
energies. The Full Width at Half Maximum is shown to be related with the 
mass concentration, suggesting that the kinetic energy of the gas in AGNs has 
a gravitational origin. Within the standard accretion model, the Transition- 
type Objects, Seyfert 2 galaxies and LINERs represent AGNs forming super- 
massive black holes 011 different mass-scales, or they could be related through 
an evolutionary process, the LINERs representing the end product of this 
evolution. 

Key Words: galaxies: activity — (galaxies:) quasars: general — galaxies: 
morphology — galaxies: accretion rate 

1. INTRODUCTION 

In the early 1960's, numerous unreso lved radio sources were found that are 
now known as quasars |Weedmanlll"986T) . The optical spectra of quasars show 
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they are extremely rich in hot gas, ionized by copious amounts of high-energ y 
photons |Osterbrocklll989t iBlandford. Netzer fc Woltierlll990l : iKrolikl 119991 ) . 



than radio-loud ones 


(Sandaee 196 


5; Braccesi et al. 


1980; 


1983; 


Padovani 


1993; 


de Vries et al. 


20061). Understanding 



jreen 



optical spectra of quasars was also instrumental in confirming the interpreta- 
tion of the redshift in terms of the expansion of the universe ( Schmidt 19631 : 
Heckman et al.lll984l ). Adopting this cosmology and putting all the obser- 



vations together, with time it became clear that the energy produced by a 
quasar is so phenomenal that the only physical mechanism known that could 
explain it is the transformation of gravitational energy into radiation through 
the accretion of m atter on a supermassive black hole (SM BH) forming at the 



center of galaxies |Reeslll978[ ISoltanl fliil : iNetzerl [l985h . Today, this inter 



pretation constitutes the basis of the standard model of the Active Galactic 
Nuclei (AGN) phenomenon. 

Subsequent observations in the 1970's and 1980's rapidly showed that 
quasars were much more common at high redshift, that is, in the past, which 



with time ( Schmic 


tlll972: Braccesi et all 1980: Veron 


1986; 


Yee & Greerll987: 


Bovle et al. 


2000: 


Richards et al. 


12006: Croom et al. 


2009). The common in- 



terpretation for this evolution is that luminous AGNs mark the first phase in 
the formation of galaxies by mergers, possibly follo wing, very tightly and in a 



complex intricate way, the formation of their bulges ( Cavaliere & Szalav 


1986; 


Silk & Rees 1 


)98: Miller et al. 2003; Di Matteo et al. 


2005: 


Laoi et al. 


2006; 


Hopkins et al. 


2007at 


Horst & Duschll 2007t 


Elbaz et al. 20091: Letawe et al. 


20ld Treister et al.l 20101). Consistent with this interpretation it was found 



velocity dispersion of the galaxies (Magorrian et al.lll998; Ferrarese & Merritt 


2000 


Gebhardt et al.l2000; 


Haring & Rixll2004: Peterson et al.l2005: Giiltekin et al. 


2009 


Graham et al. 20 111). 



So far, the standard model for AGN seems to offer an elegant and credible 
explanation for quasars, in good agreement with cosmology. However, here 
lies one mystery: if one counts the number of luminous AGNs at high redshift, 
then almost all massive galaxies in the nearby univer se must have had a com- 
parable high-activity phase during their formation ( Kormendv fc Richstonel 



1995; Richstone et al 



1998) 



But what is the evidence at low redshift in fa- 
vor of such a conclusion? In other words, where are the luminou s quasars, 
or rather their remnants, at low redshift (|L6pez-Corredoiral l201lh ? Within 
this interpretation, the luminous AGN phase is assumed to be very short in 
duration (10 6 to 10 8 yrs) compar ed to the formation and evolution of galaxies 
( Martini 2004 : Kelly et al. 201dh . This suggests that most quasar remnants 
today should not show any tr ace of non-stella r activity in their nucleus. But 
how short really is this phase ( Cattaneol 2001 ) , and how consistent is this sce- 
nario wit h the low level of nuclear activity observed in galaxy surveys at low 
redshift JlJeda et al.l 120031: iMerlonil liool iHeckman et~al1l2004 iBarger et af 
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2005))? 

In the present study we use one of t he largest spectroscopic survey to 
date the Sloan Digital Sky Survey (SDSS; lYork et alJbood : Istoughton et al 



20021 ) , in order to shed some light on the evolution of AGNs at low luminosity. 



The SDSS spectral survey has revealed that a high fraction of nearby galaxies 
possess ionized gas in their nucleus. A small fraction of these galaxies, the 
Seyfert 1 (S yl), show broad emission- line components akin to what is ob served 



in quasars |weedmanlll986l |Pibail[l987i: IOsterbrocklll989t lKroliklll999h . The 



standard model of AGN suggests that the S yls are similar to quasars, formin g 
the low-luminosity tail of their distribution (jSmith et al 1 ll98dlMarshalill98i . 
In terms of the SMBH theory it is expected that these AGNs have lower-mass 
black holes for the sam e accretion luminosity , L acc , and efficiency, 77, than the 
more luminous AGNs ([Hopkins et al.ll2006al ). 

However, much more numerous in the SDSS survey are narrow emission- 
line galaxies (N ELGs) . Using diagnostic diagrams that compare va rious emis- 



sion-line ratios (Baldwin et al 



1981 



rent classification criteria (IKewlev et al 



Veilleux fc Osterbrocklll987l). a nd diffe- 
Kauffmann et al 



2001 



fraction of these NELGs were found to have spectra consistent with that of 
an AGN, although they do not show as broad lines as the Syl and are less 
luminous. One of these AGNs, the Seyfert 2 (Sy2), is thought to differ from 
the Syls only because the broad-line regions in these galaxies are obscured 
by a ci rcumnuclear torus of gas and dust (this is the so-called unified AGN 
model; lAntonucci fc Millerlll985l) . On the other hand, the Sy2s may also be 
intrinsically different from the Syls, lacking, for some unknown reason, the 
broad-line regions. 

Another type of AGN which is even more proble matic than Sy2 is the 
Low Ionization Nuclear Emission-line Regions (LINER; lHeckmanlll980l : ICozioll 
19961: IKewlev et al.ll2006l) . Despite the classification based on standard spec- 



tral diagnostic diagrams, which identify these galaxies as AGNs similar to 
Sy2s, many researchers today claim that LINERs are "false AGNs", mostly 
because many show active and intense star-formation activity in their circum- 
nuclear regions, and do not present clear evidence for a compact source in 
their nucleu s. The presence of star f ormation in narrow-line AGNs seems now 
indubitable (jKauffmann et al.ll2003l) . This is why the Transition- type Object 
(TO) class was systematically included, defining in diganostic diagrams a sort 
of "buffer" zone between pure Star Forming Galaxies (SFGs) and pure (or 
rather dominant) AGNs. The presumed absence of a compact source in LIN- 
ERs, on the other hand, was never proven. Instead, evidence for an AGN 
contribution is usually downgraded in favor o f star formation, but never elim- 
inated completely (e.g.. lAnnibali et al.ll2010h . 

Finding evidence of SMBHs ge nuine AGNs s hould not pose such a problem. 
For example, in one of his books, Osterbrock ( 19891 ) proposed a simple test, 
which consists in comparing the luminosity of the emission line, say Ha, with 
the continuum luminosity in the blue, say at 4800A (sufficiently far from 
H/3 at 4861 A). According to the standard model of AGN, the source of the 
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ionizing flux is the UV continuum produced by the accretion of matter onto 
a SMBH. Consequently, a typical AGN must show a characteristic power law 
relation, L em oc L" ont , significantly different from any relation followed by 
SFGs (the same principle, for the same reasons, is at the basis of the spectral 
diagnostic diagrams). Up to now, however, this test was never applied to a 
large and statistically significant sample of NELGs (e.g., Buson et al. 1993t 



Macchetto et al.lll996f ). This is one of the major goals of our study. 



The plan of this article is the following. In Section 2, we describe the selec- 
tion of our sample of NELGs and explain how the data for our analysis were 
obtained. In Section 3 we present the results of our spectral classification anal- 
ysis. In Section 4 we study and discuss the results of the Osterbrock power-law 
test and in Section 5 we present other evidence consistent with SMBHs at the 
center of NELGs. Finally, in Section 6 we discuss how the standard accretion 
model can explain the different types of AGNs. Our conclusions can be found 
in Section 7. 



2. DESCRIPTION OF THE SAMPLES AND DATA FOR ANALYSIS 
The data for our study come from the m ain catalogue of the Sloan Digita l 



Sky Survey, Data Release 7 (SDSS, DR 7; lAdelman-McCarthv et all 120071) . 
Using the Virtual Observatory (VO) servic^f), we downloaded the spectra 
of 229618 galaxies, all classified as NELGs, and with the highest spectral 
quality, namely emission lines with signal-to-noise ratio S/N > 3 on H/3, 
[OIII]A5007A, [NII]A6584A and Ha and S/N > 10 for the adjacent continuum. 
To minimize the Malmquist bias effect (and avoid problems with cosmology), 
we have limited our sample to a redshift z < 0.25. 

In the VO, the spectra are corrected for Galactic extinction, shifted to their 
respective rest frame, resampled to AA = lA between 3400 and 8900A, and 



proce ssed by the spectral synthesis code STARLIGHT (jCid Fernandes et al 



20051 ). The STARLIGHT code produces automatically a stellar template from 
which important characteristics of the NELG hosts can be deduced. In par- 
ticular for our study, we retrieved from the templates the velocity dispersion 
of the stars, a ap , in the bulge of the galaxies. It is important to mention 
that STARLIGHT does not allow to correct for rotation. However, this prob- 
lem is reduced in our study because the small aperture of the SDSS fiber (3 
arcseconds) projected at the distance of the gala xies always covers a region 



well w ithin the bulge, which is pressure supported (jBender. Burstein fc Faber 



19921 ). We have also verified that within the resolution of the stellar libraries 
used in our study a ap does not change significantly as the aperture increases. 
The velocity dispersion estimated by STARLIGHT converges rapidly to the 
real stel lar velocity dispersion of the bulge, cr*. This is consistent with the 



study of iBernardi et all (|2003al) , who have shown that the measured using 
SDSS spectra have an uncertainty smaller than 14%, which is equal to the 
intrinsic variance of the velocity dispersion as measured in galaxies having 
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the same morphology and redshift. Therefore, using the effective radius, R ei 
and a ap the mass of the bulge of the galaxy in our sample can be reasonably 
estimated by applying the virial theorem. Once the bulge mass is known, we 
can i nfer the mass of th e SMBH at the center of the AGNs in our sample 
(e.g.. lHaring fc Rixll2004l) . Howev er, for our study we decided to use instea d 
the relation between Mbh and ( Gultekin et al. 2009t Graham et al. 2011 ). 
which is more direct because it does not require determining the effective radii 
of the galaxies. 

To separate the NELGs in our sample according to their different nuclear 
activity, we used three different spectral-line diagnostic diagrams. The line 
ratios necessary for the classification were measured in the corrected spectra, 
obtained by subtracting the stellar templates produced by STARLIGHT from 
the raw spectra. This method improves greatly the classification process by 
reducing the diluting effect on emission-lines of the underlying older stellar 
populations. On the other hand, some uncertainties are introduced by this 
method because STARLIGHT assumes the continuum is purely stellar. How- 
ever, this affects more the characteristics of the stellar templates than the 
emission line ratios. That is, in AGNs the velocity dispersion and the ages 
of the stellar population may be slightly underestimated. Note however that 
to apply the Ostcrbrock test, the luminosity of Ha is determined using the 
corrected spectra, which is free of stellar contribution, while the luminosity 
of the continuum, at any wavelength (we used two in our study) , is measured 
in the raw spectra, which include the contribution of stars and the AGN if 
present. 

To verify if the power laws in our test depend on the morphologies, we 
have also determined this parameter for all the gal axies in our samp l e. To 
do so, we have followed the method developed by IShimasaku et al.l (|200ll ) 



and lFukugita et al.l (|2007l ). Using the photometric colors, u — g, g — r, r — i, 
and i — z, as defined in the ugriz photometric system of SDSS^, and the 
inverse conc entration indice s, i?5o(r)/i?g (r), which is the ratio of the Pet- 
rosian radii ( Petrosian 19761 ) containing 50% and 90% of the total flux in the 
r band, c orrelations were established with the s tandard Hubble morpholo gi- 
cal types (|de Vaucouleurs et al.lll99ll ). Following iBlanton fc Roweisl (|2007l ). a 
K-correction was applied on the magnitudes. 



3. NUCLEAR ACTIVITY CLASSIFICATION 



As we have mentioned in the previous section, our method allows us to 
obtain a very precise classification by activity types. In Figure [T] we present 
our primary classification based on the most useful of all of the diagnostic 
diagr ams (the BPT-VO diagram; iBaldwin et alJll981[ IVeilleux fc Osterbrock 
19871 ). The separation crite ria applied to disting uish between AGNs, TPs an d 
SFGs are those proposed bv lKewlev et al.l (|200ll ) and lKauffmann et al.l (|2003f ). 
Results for our classification are shown in Table Q] There we can see that the 



http:/ /casj obs.sdss.org 



G 



TORRES-PAPAQUI ET AL. 



SFG type is the dominant activity in NELGs, with 53.5% of the sample. The 
remaining 46.5% are separated into 28.8% TOs and 17.7% AGNs. As we 
mentioned earlier, these results are consistent with the view that AGNs are 
rare at low redshift, and frequently accompanied by star formation activity 
in the circumnuclear regions, which greatly complicates their identification 
and/or study. 

| | i i i i | i i i I | i i i J | i i i r— | i—i I 




AGNs 



-1.5 -1.0 -0.5 0.0 0.5 

Log ([Nll]X6584/Ha) 



Fig. 1. BPT- VO diagram for th e NEL Gs. The criterion that separates SFGs from 



TOs is that of 
TOs is that of 
and the AGNs in red. 



Kauffmann et al 



( 20031 ) and the criterion that separates AGNs from 



Kewlev et alJ ( 2001 ). The SFGs are shown in blue, the TOs in green 



It is important to emphasize that despite having determined the line ratios 
with great accuracy, the BPT-VO diagram still shows a continuous distribu- 
tion. Therefore, the separation criteria applied in this diagram look some- 
what arbitrary, allowing numerous crossed activity identifications or ambigu- 
ous classifications. This is another important difficulty in studying narrow-line 
AGNs, which could have serious consequences, in particular when the samples 
studied are small. 
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TABLE 1 

NUMBER OF GALAXIES IN SAMPLES 





SFGs 


TOs 


AGNs 


Total 


122751 


66234 


40633 


SFGls 


67513 






SFG2s 


55238 






[01] Sy2s 






12264 


[01] LINERs 






14102 


[SII] Sy2s 






12911 


[SII] LINERs 






19893 


Sy2s 






16970 


LINERs 






23663 



The difficulty becomes more obvious when one tries to differentiate be- 
tween LINERs and Sy2s. No natural distinction appear s in the BPT-VO di- 
agram between these two activity types. In their study, Kewlev et al. ( 20061 ) 
suggested that LINERs and Sy2s may be distinguished based on differences in 
the [OIJA6300 and [SII] AA6717, 6731 emission-line intensities. This is shown 
in Figure [5] where we reproduced the analysis of these authors, applying their 
proposed criteria on our sample of AGNs to separate Sy2s from LINERs. The 
results are reported in Tabled] Note that only 65% of the AGNs in our sample 
can be classified using [OI]A6300, which is a consequence of the intrinsically 
low intensity of this line in these galaxies (affecting LINERs more than Sy2s). 
On the other hand, even when using a lower-ionization line like [SII] only 
81% of the AGNs can be classified this way. Moreover, some galaxies change 
activity type passing from one diagram to another. 

To be able to classify all the AGNs in our sample, we used the results 
of the analysis in Figure [2] to establish a dividing line between the LINERs 
and Sy2s in the BPT-VO diagram. This is shown in Figure [3] In the right 
panel in b) we have located in the BPT-VO the galaxies classified as Sy2s in 
Figure [2] In d) we have located those classified as LINERs, and in c) we have 
traced the position of the galaxies that change activity type depending on the 
diagnostic diagram used. In a) we added the three subsamples together, which 
represent 53% of the whole AGN sample. Based on this result, it is clear that 
we cannot eliminate completely the arbitrariness in the separation between 
the LINERs and Sy2s in any diagnostic diagram. As a compromise, we chose 
a dividing line that splits the sample of ambiguous types in c) in two groups 
having the same number of galaxies and used this separation as our principal 
classification criterion in the BPT-VO diagram to distinguish between Sy2s 
and LINERs. The separation corresponds to the following linear relation in 
Log: 
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Fig. 2. Diagnostic diagrams distinguishing between Sy2s and LINERs (jKewlev et al 
200d ): a) [OIII]/H/3 versus [SII]/Hq; b) [OIII]/H/3 versus [OI]/Ha. 



Log 



([oni}\ 
V hp ) 



= 1.64 x Log 



m 



0.38 



(1) 



Note that our chosen criterion slightly favors Sy2s over LINERs (that is, 
in Figure [31 a slightly higher number of LINERs in d) than Sy2s in b) change 
activity classification). After applying this criterion to all the AGNs in our 
sample we count 58% LINERs and 42% Sy2s (see Table HJ. 

The presence of star formation in narrow-line AGNs is now well recog- 
nized. However, this was not always the case in the past. At the time when 
the diagnostic diagrams were developed, the main assumption was that the 
right part of the continuous ^-shaped distribution traced by the line ratios in 
the BPT-VO diagram was occupied by pure AGNs. With the introduction of 
the TO type the preoccupation now is slightly different. The current question 
is whether there could be AGNs among the SFGs. For instance, in Figure Q] 
some SFGs show unusually high emiss ion-line ratio s of [Nil] /Ha as compared 
to what is observed in HII regions (see lCozioIl996h . This can be explained ei 



ther by assuming an excess of excitation due to an AGN (e.g.. lStasinska et al 



20061 ) , or by assuming an excess of nitrogen abu ndance, related with the for 



mation of massive bulges (|Coziol et al.lll999h . In lTorres-Papaqui et al.l (|2012T ) 



SMBH IN NELGS 



9 




-1.0 -0.5 0.0 0.5 -i.o -0.5 o.o 0.5 

Log ([Nll]X6584/Ha) Lo 9 ([Nn]X6584/Ha) 



Fig. 3. BPT-VO diagram for the NELGs classified as Sy2s and LINERs in Figure [2] 
in b) we show all the galaxies classified as Sy2s using two diagnostic diagrams; in 
d) we show those classified as LINERs; and in c) we show those that change type 
according to the diagnostic diagram. The continuous line separates this last sample 
in two equal-sized subsamples. This dividing line is reported in b) and d) and on 
the whole sample in a). 



this problem was explored in great detail and it was concluded that, in gen- 
eral (except for a few ambiguous cases, as expected based on the ambiguity 
of selection criteria due to the continuous distribution of the NELGs in the 
BPT-VO diagram), there is no evid ence of AGNs in the SFG s, in good agree- 
ment with the criterion proposed by KaufFma nn et al.l (|2003l) to separate TOs 
from SFGs. Following iTorres-Papaaui et al.l (|2012l ). we have separated our 
SFG sample in two. The SFG1 (55.2% of the SFGs), are consistent with the 
standard view of late-type spirals having experienced constant star formation 
over their whole history. The SFG2s with an excess of nitrogen abundance 
form 44.8% of the SFGs. These are massive early-type spiral galaxies (that 
is, spirals with massive bulges) that have experienced (or some are still expe- 
riencing) intense bursts of star formation during their formation. 

In Figure 2] we show the final separations in activity type in the BPT- 
VO diagram, according to which all the galaxies received only one classifi- 
cation. The number of galaxies in each sample is identified in bold in Ta- 
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Fig. 4. BPT-VO diagram defining the five different activity types used in our study: 
SFG 1 (blue), SFG 2 (black), TO (green), Sy2 (red) and LINER (magenta). 



blc [T] There are almost as many TOs as there are SFGls. The number 
of SFG2s is also remarkably high, which indicates that powerful starburst 
galaxies are not rare in the nearby universe. The "pure" AGNs on the 
other hand are much less common. Among these, the LINERs arc domi- 
nating over the Sy2s (this is despite the bias introduced by our classification 
criterion), which is consistent with previous surveys, showi ng that LINERs 
are very comm on in intermediate spiral galaxies in the field ( Heckman 1980t 



Ho et al.lll997l ). and constitute the pri ncipal activity type of galaxies encoun 



tered in dense galactic environments (Phillips et al. 1986 : Coziol et al. 19981 : 



Miller et al 



200S 



Martinez et al.l 120081 12010 1 Note that by extending the 



separation between LINERs and Sy2s to the TO region, more than 2/3 would 
classify as TO/LINER rather than TO/Sy2 galaxies. This result reflects the 
ambiguity found in the literature about the real nature of LINERs. 



SMBH IN NELGS 



11 




27 28 29 30 27 28 29 30 
Log (L re U800 / ers s -1 Hz -1 ) Log (L FC X4800 / ers s" 1 Hz -1 ) 



Fig. 5. Results of the Osterbrock test for the NELGs separated by activity type: 
top-left, the SFGls are plotted over the SFG2s (almost no difference in distribution); 
top-right, TOs; bottom-left, Sy2s and bottom-right, LINERs. The parameters of the 
linear fits are given in Table [2] 



TABLE 2 

PARAMETERS OF LINEAR FITS AND CORRELATION 
COEFFICIENTS 



Sample 




P 




r 


SFGs 


1.09 ± 0.04 


9.90 ± 0.35 


0.91 


0.91 


SFGls 


1.10 ± 0.05 


9.66 ± 0.34 


0.88 


0.89 


SFG2s 


1.14 ± 0.03 


8.41 ± 0.34 


0.93 


0.93 


TOs 


1.16 ± 0.04 


7.51 ± 0.37 


0.77 


0.79 


Sy2s 


1.48 ± 0.03 


-1.31 ± 0.03 


0.81 


0.82 


LINERs 


1.45 ± 0.03 


-1.21 ± 0.03 


0.75 


0.75 
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4. THE OSTERBROCK TEST: FIRST EVIDENCE FOR THE 
PRESENCE OF A SMBH IN NELGS 

The Osterbrock test is based on the principle that the source of ionization 
of an AGN is the continuum produced by the accretion of matter on the 
SMBH at its center. By comparing the luminosity of the line-emission with the 
luminosity of the continuum in the blue part of the spectrum, the AGNs should 
therefore follow a characteristic power law that can be easily distinguished 
from any relation produced by other ionizing sources, in particular O and 
B stars like in the SFGs (and possibly TOs). For our test we compare the 
luminosity in Ha with the luminosity of the continuum at 4800A. The results 
are shown in Figure [5] for the NELGs in our sample separated by activity type. 
In Table [2] we give the parameters of the linear fits and the Spearman rank 
(r s ) and Pearson (r) coefficients of correlations. All the NELGs show strong 
linear correlations in Log of the kind: 

T /L(Ha)\ T /L FC (4800A)\ n 

Log = a * Log tCA _ + /3 2 

\ erg s 1 / \ erg s 1 Hz 1 J 

Examining the parameters of the linear fits in Table [2j we see that the 
SFGs and TOs have similar a and /3, suggesting their dominant source of 
ionization is similar. The LINERs and Sy2s also have similar a and j3, which 
are significantly different from those determined for the SFGs. The trend 
is for the slope to be steeper in AGNs than in SFGs. That is, for the same 
luminosity in the emission line, the luminosity of the continuum is higher in an 
AGN than in a SFG. This is the main characteristic of AGNs. Therefore, we 
conclude that the ionization source in all the LINERs in our sam ple is the same 



as in the Sy2s. Consequently, both types are genuine AGNs (jKewlev et al 



20061 ). 

Comparing the results of the linear fits in Figure [6] a more complete and 
consistent picture can be drawn. The LINERs are not only similar to the 
Sy2s by their ionization source, but also have significantly lower emission-line 
luminosity. In other words, the LINERs are the low-luminosity equivalent 
of the Sy2s, or better, since the LINERs are much more common than the 
Sy2s, the Sy2s are the high-luminosity equivalent of the LINERs. At the 
same time, we note that the TOs, although similar to the SFGs, present the 
same tendency of a steepening of the slope, which seems consistent with the 
standard interpretation of a TO as a galaxy having a double nature, that is, 
the TO type is a mixture of SFG and AGN. 

The Osterbrock test is unambiguously suggesting that the ionization source 
in the narrow-line AGNs is the continuum produced by the accretion of mat- 
ter onto a SMBH at their center. Therefore, the power law should be the 
same, independent of the morphology. In Figure [7] we show the results of 
our morphological classification for all the galaxies in our sample separated 
by activity types. There it can be seen that the morphologies change signifi- 
cantly from one activity type to the other. The SFG2s, for example, reside in 
spiral galaxies with slightly earlier type (Sbc) than the SFG1 (Sc/Scd). The 
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Fig. 6. Comparison of the linear fits for the luminosity power law in NELGs having 
different activity types. 

TOs (Sab/Sb/Sbc) have even earlier type than the SFGs. They are, on the 
other hand, almost indistinguishable from the Sy2s. The LINERs have the 
earliest types of all (SO/SOa). 

To check if the linear relations varied with the morphology of the galaxies 
we have separated the morphology distributions in three different groups. The 
results of the linear fits in these groups arc presented in Tabic |3l As expected, 
we observe no significant variation of the slope with the morphology. This 
suggests that the variations in stellar populations, being significantly older in 
the early-type galaxies, do not influence the linear relation between the ionized 
flux and continuum. This shows that the slope of the linear relation depends 
only on the activity type. In the SFGs and TOs young massive stars, which 
are common in all star-forming galaxies, are dominating the continuum in the 
blue part of the spectrum, explaining the similar slopes. Similarly in the Sy2s 
and LINERs, the same characteristic slope suggests a common phenomenon 
different from star formation, which is the accretion of matter onto a SMBH 
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TABLE 3 

PARAMETERS OF LINEAR FITS AND CORRELATION 
COEFFICIENTS ACCORDING TO ACTIVITY TYPE AND 
MORPHOLOGY 
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at their center. The LINERs and Sy2 are similar types of AGNs, differing 
only by their intrinsic line-emission luminosity. 



5. CONCENTRATION OF MASS IN THE CENTER OF GALAXIES AND 
FWHM: OTHER EVIDENCE IN FAVOR OF A SMBH AT THE 
CENTER OF NELGS 

The Osterbrock test for the NELGs is consistent with the presence of an 
accreting SMBH at the center of the Sy2s and LINERs. We have already 
shown that the LINERs reside in galaxies having the earliest morphological 
types of all the NELGs in our sample, which suggests they have generally 
massive bulges. However, the TOs show no important morphological differen- 
ces from the Sy2s, which suggests that the presence of a massive bulge is not 
sufficient to explain why a SMBH becomes active in a galaxy. A more restric- 
tive constr aint is the gravitational binding energy (or gravitational potentia l 
energy; see 



Hopkins et Zll2007bt Uller fc Richstonell2007l; ICoziol et allboill) 



Assuming the formation of a SMBH follows the formation process of galaxies 
by gravitational collapse, then galaxies developing such objects at their nu- 
cleus should have naturally higher binding energy than those without it, and 
this would be irrespective of the morphology of the galaxy. By definition, a 
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TABLE 4 

VMA AND FWHM IN NELGS SEPARATED ACCORDING TO 
ACTIVITY TYPE 



Log VMA Log FWHM(Ha) 
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Fig. 8. Box-whisker plots for VMA and FWHM. The extremely small notches (barely 
visible) are due to the large sizes of the samples, and reflect the significant differences 
in medians of these samples. The dots correspond to the means. 



higher gravitational binding energy implies that more mass is being concen- 
trated at the center of these galaxies. Having at hand the velocity dispersion 
in the aperture of the spectra, a ap , we can easily estimate the virial mass 
within the aperture (VMA), and use this as a proxy for the binding energy: 
the higher the mass within the aperture the higher the binding energy. 
We define the VMA as: 

VMA = ^ (3) 



where R is the physical radius projected by the aperture. The box-whisker 
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plots in Figure[8£i reveal significant differences of the VMA as a function of the 
activity type. The values describing the statistical distributions are reported 
in Table |U The median and mean for the LINERs and Sy2s are significantly 
higher than for the SFGs and TOs. An AN OVA statistical test confirms the 
significance of these differences at a 95% confidence level (sec appendix A). 
From our analysis we find that the VMA gradually increases from the SFGls 
and SFG2s, to the TOs and Sy2s, culminating in the LINERs. 

In Figure[9]we show the median VMA for the NELGs with different activity 
types and different morphologies. The VMA is always higher in the AGNs. 
The fact that this characteristic is independent of the morphology of the 
galaxies is important. It clearly shows that the VMA does not depend on the 
mass of the bulge, as we suggested. The AGNs have higher VMA because 
they have higher gravitational binding energies, which is due to the formation 
of a SMBH at their centers. 

Note also the different trends in Figure [9] the LINERs have significantly 
higher VMA than the Sy2s in early-type galaxies, but the difference between 
the two decreases in the later types, while the Sy2s differ more significantly 
from the TOs in the late types than in the early types. In general, the TOs 
have a VMA > 10 10 M Q , significantly higher than the SFG2s. Since they have 
physical characteristics comparable to those of the Sy2s, differing possibly only 
in their higher star-formation rates, it seems plausible to assume that a SMBH 
is also present at the center of these galaxies. 




Fig. 9. Median VMA as a function of morphology, for galaxies with different activity 
types. 
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According to its definition (cf. Eq. 3), the VMA varies linearly with the 
projected aperture size, and consequently it varies with the redshift. We show 
in Figure [10] the median VMA as a function of redshift for galaxies having 
different activity types. It can be seen that at any redshift the VMA is always 
higher in the AGNs. Note that we do not observe much difference between the 
TOs and Sy2s in this diagram, which reinforces the idea that TOs also have 
a SMBH at their center. From this comparison it is clear that AGNs have a 
characteristically high VMA, independent of the the projected aperture size. 
That is, going from high to low redshift the projected aperture size decreases, 
but the VMA of an AGN keeps being higher than in other activity types. At 
the limit where the projected aperture size become unresolved, this behavior 
is consistent with the presence of an unresolved massive object at the center 
of all the AGNs. 



1 1.0 




0.05 0.10 0.15 0.20 0.25 

Redshift (z) 

Fig. 10. Median VMA as a function of redshift for galaxies with different activity 
types. 



5.1. Relation of VMA with the FWHM 



In the past (e.g., lOsterbrockl 119891 ) . it was fairly well accepted that one 
defining characteristic of AGNs is that their emission-lines have large Full 
Width at Half Maximum (FWHM). Narrow-line AGNs are no exception. For 
instance, Sy2s have permitted and forbidden lines of similar widths in the 
range 500 km s , which is significantly broader than in a typical SFG where 
the FWHM is of the order of 200 km s~ x . Here we show that this characteristic 
is directly related with the VMA, and a consequence for AGNs of having an 
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Fig. 11. Relation of FWHM with VMA. The points correspond to the medians for 
different activity types and the bars cover the first and third quartiles. 

unresolved SMBH at their center. In FigureUJa, we show the box-whisker plots 
for the FWHM of the Ha line in galaxies having different activity types. We 
observe that the FWHM increases in the same way as the VMA. This can also 
be observed in Table |H As for the VMA, an AN OVA statistical test confirms 
that the differences are statistically significant (see appendix A). In Figure QT] 
we report the VMA and FWHM medians, first and third quartiles for the 
NELGs with different activity types. The dashed line is a linear relation that 
quantifies the gradual increase of VMA with the FWHM: 

Log (™) . (3.39 ± 0.03) x Log (™?Si) + d-4 ± 04) (4) 

Based on this result, we see that the TOs with FWHM > 10 2 5 km s _1 are 
more similar to Sy2s than to SFGs. 

The reason why AGNs have characteristically large FWHM is not well 
understood. In terms of dynamical Doppler effects, broad-line profiles trans- 
late into high-speed gas motions, either rotational and/or turbulent. What 
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is clear, on the other hand, is that the high level of energy necessary to pro- 
duce these movements cannot co me from thermal sources, namely massive 
starburst winds or SN explosions (jBlandford. Netzer fc Woltierlll990n . which 
anyway seems difficult to assume in LINERs where the star formation activity 
may be at its lowest level. Our result suggests that this energy is gravitational, 
related with the high binding energy produced by the formation of a SMBH 
at the center of the galaxies. 



6. DISCUSSION: MASSES OF THE SMBHS AND ACCRETION RATES 

IN NARROW-LINE AGNS 
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Fig. 12. Box-whisker plots for the masses of the SMBHs and accretion rates in the 
narrow-line AGNs. The notches, which are barely visible due to the large sizes of 
the two samples, reflect the significant differences in medians between the samples. 
The dots represent the means. 

In the two previous sections we presented new evidence for the presence 
of SMBHs at the center of NELGs. Added to those already recognized in 
the literature (in particular the BPT-VO diagnostic diagram) this evidence 
strongly supports the interpretation that narrow-line AGNs, like the Sy2s 
and LINERs, are genuine. Within this interpretation, we now discuss how 
the standard AGN model, the accretion of matter onto a SMBH, may explain 
the differences between the different types. For this discus sion, the masses of 
the SMBHs were estimated following Giiltekin et al. ( 20091) and the accretion 



rates correspond to l^Boi/^Edd, where the bolometric luminosity is obtained 
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TABLE 5 

SMBH MASSES AND ACCRETION RATES IN NARROW-LINE AGNS 
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Median Mean 
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Sy2s 


7.25 


7.54 


7.53 


7.80 


-1.50 


-1.29 -1.28 


-1.06 


LINERs 


7.55 


7.86 


7.87 


8.20 


-1.92 


-1.65 -1.65 


-1.41 



using the approximate relation Lb i = 9ALa(5100A) (jKaspi et al.ll2000h . The 
continuum is evaluated at a different wavelength than before in order to be 
able to co mpare our results with th ose of Broad-Line AGNs (BLAGNs; in 
particular beterson et alj|2004 120051) . 

In Figure [12] we draw the box- whisker plots for the masses of the SMBHs 
and their accretion rates. The notches once again are barely visible due to 
the large sizes of our samples, also emphasizing the statistically significant 
difference of the medians. The values describing the distributions are reported 
in Table [SJ The statistical significance for the differences of the means was 
established using ANOVA t ests (see appendix A ). Our results are comparable 
to those obtained before bv lKewlev et al.l (|2006l ). who used different methods 
to obtain the masses of the SMBHs and accretion rates. Our main conclusion 
is consequently the same as these authors: the LINERs have more massive 
SMBHs that accrete at lower rates than the Sy2. 

In Figure [T31 we compare the masses of the SMBHs and the accretion 
luminosity as we measure d in the Sy2s with th ose estimated i n BLAGNs ob 



served in reverberation bv lPeterson et al.l (|2004l) and studied in lPeterson et al 



(|2005l ). In general the Sy2s have lower- mass SMBHs than the BLAGNs. In Fi- 
gure [13] the Sy2s seem also to be accreting at lower rates than the BLAGNs, 
i.e., below 0.1. Note that for comparable masses of SMBH the accretion 
rates we have determ ined are in good agreement with those estimated by 



Peterson et al.l (|2004l) . Our much larger sample however allows us to observe 



a new trend for the more massive SMBHs in the Sy2s to have lower accretion 
rates. If we trace a linear correlation over the data the slope for the Sy2s 
is above 1, 1.17 ± 0.04 to be exact, with correlation coefficients of 0.6 (both 
Spearman rank and Pearson, with probabilities of obtaining the same val- 
ues by chance practically zero). This correlation suggests that as the SMBH 
increases in mass, its luminosity, or accretion rate, decreases. 

In Figure Owe compare the LINERs with the BLAGNs. Here the diffe- 
rences are much more important. In the LINERs some of the SMBHs reach 
almost 10 9 Mo, comparable to the most luminous BLAGNs. However, they 
accrete at much lower rates, around 0.01. Like for the Sy2, but with a steeper 
slope of 1.43 ± 0.04, a linear correlation (also with correlation coefficient 0.6) 
suggests that as the SMBH increases in mass the accretion rate decreases. 

In parallel with the variation of accretion rate with the SMBH mass, we 
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Fig. 13. Mass of SMBH vs. accretion luminosity for the Sy2s. We show the data as 
shaded density contours; the darker the region the more numerous the data points. 
The indivi dual data points cor respond to the sample of BLAGNs as previously 
studied bv lPeterson et al. I (|2005l ). The continuous lines are different accretion rates, 
Lsoi/Lfidd- The dashed line is a linear fit to our data. 



also observe a variation in morphology. In Figure [T5] for the Sy2s we show 
the SMBH mass as a function of accretion luminosity in three different mor- 
phological groups: in a) the early-type group Tl, in b) the intermediate- type 
group T2, and in c) the late- type group T3. The definition of these morpho- 
logical groups are given in Table [3] We see systematic changes, the SMBH 
masses increasing from the late types to the early types, and the accretion 
rates decreasing in the same sense. Moreover, the slope of the linear corre- 
lation increases from 1.09 ± 0.04 in the group T3 to 1.13 ± 0.03 in the T2 
and Tl groups. The same trend is observed for the LINERs separated in 
morphological groups in Figure [TBI There is a clear augmentation in the mass 
of the SMBH passing from the late types to the early types, concurrent with a 
significant decrease in the accretion rates. The slope of the linear correlation 




Fig. 14. Same as Figure H51 for the LINERs. 



also steepens passing from 1.21 ± 0.03 in the T3 to 1.34 ± 0.05 in the T2 and 
1.31 ± 0.04 in the Tl. 

In Table |6] we give the characterictic statistics for the masses of the SMBH 
as they varied with the morphoplogcal types and in Table [7] we give the char- 
acterictic statistics for the luminosities. For the Sy2s, on average a small 
dimming in luminosity by 18% is observed going from the T3 to the Tl, but 
not in the LINERs where the luminosity increases by 11%. The differences 
in luminosity between the LINERs and Sy2s are also marginal: the Tl LIN- 
ERs are 15% more luminous than the Tl Sy2s, while the T3 are 13% dimmer. 
However, at the same time we see a constant increase in the SMBH mass from 
the T3 to the Tl in both the Sy2s (the Tl are 53% more massive than the 
T3) and the LINERs (the Tl arc 75% more massive than the T3), and we see 
an increase by 70-85% comparing the LINERs with the Sy2s by morpholo- 
gies. Therefore, for almost the same bolometric luminosity the mass of the 
SMBH increases significantly with the morphological type, which translates 
into a decrease in accretion rate (or accretion efficiency) with the mass of the 
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Fig. 15. Same as Figure [13] separating the sample of Sy2s in three different mor- 
phology groups: a) early-type Tl; b) intermediate-type T2; c) late-type T3. The 
definition of the groups can be found in Table [3] 
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Fig. 16. Same as Figure IT51 for the LINERs. 



SMBH. 

The above trends are consistent with some pattern in evolution (e.g. lKellv et al 



20101 ) . As the bulge of the galaxies grows (due to star formation) the SMBH 



grows in mass through accretion. The bulges stop growing when no more stars 
form, which means the amount of available gas in these galaxies is already too 
low, and possibly also the amount of gas to be accreted by the SMBH, explain- 
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TABLE 6 

MASSES OF SMBHS IN NARROW-LINE AGNS HAVING DIFFERENT 
MORPHOLOGICAL TYPES 
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Log (M BH 
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Tl 


7.44 


7.67 
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7.38 


7.38 
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TABLE 7 

LUMINOSITY OF NARROW-LINE AGNS HAVING DIFFERENT 
MORPHOLOGICAL TYPES 
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Log 
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Tl 


43.16 


43.39 
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43.58 
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43.20 
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43.63 




T3 


43.21 


43.45 


43.43 
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Tl 


43.25 


43.43 


43.42 


43.60 


LINERs 


T2 


43.20 


43.41 


43.39 


43.60 




T3 


43.17 


43.38 


43.37 


43.56 



ing the lower accretion luminosity. Note that for the Sy2s with accretion rates 
near 0.1 liEdd in the group T3, the SMBH may continue to grow significantly 
in mass, and if this accretion is accompanied by intense star formation, the 
bulge would also grow, eventually transforming the morphology of the galaxy 
to an earlier type. Consistent with the relation between the SMBH mass and 
the mass of the bulge, our observations suggest the Sy2s are still very much 
active in building both, their bulges through star formation concurrently with 
their SMBH throuh accretion. 

Different from the Sy2s, the accretion rates in the LINERs are much lower, 
which suggests that the LINERs may already be at, or very near, the end 
of their formation process, which is also consistent with their earlier mor- 
phologies. Therefore, the LINERs may have been more active in the past. 
Considering the similarity in binding energy and variation with morphology 
(Figure [9]) this higher luminosity phase for the LINERs may have been similar 
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to that in the Sy2s. Consequently, it may be that in each group of morphology 
the Sy2s will eventually evolve with time into something similar to a LINER. 

7. CONCLUSIONS 

We have shown that the NELGs, identified as LINERs and Sy2s, trace 
similar power laws between the luminosity in emission and the luminosity 
in the continuum, the LINERs differing from the Sy2s only by their lower 
luminosity. This result is consistent with the standard model of AGN, which 
states that the principal source of ionizing photons in AGNs is the accretion 
of matter onto a SMBH at their center. Also consistent with the presence of a 
SMBH at the center of NELGs, we have shown that the VMA, that is, the mass 
concentrated at the center of the galaxies, is always higher in the AGNs than in 
galaxies showing a higher level of star formation, like the TOs and SFGs. The 
fact that this is independent of the morphology of the galaxies suggests this 
characteristic is not merely related to the formation of the bulge, but to a more 
fundamental aspect of the formation of the galaxies. This higher concentration 
of mass in the AGNs is related with a higher gravitational binding energy due 
to the formation of a compact massive object, a SMBH, at their center. We 
have also found a strong correlation between the VMA and the FWHM, which 
suggests that the energy necessary to produce the large movement of gas in 
the AGNs is gravitational in nature, related with the binding energy of the 
SMBH at their center. 

We conclude that the standard AGN model depicts in an extremely con- 
sistent way the narrow-line AGN phenomenon observed at low redshift. Ac- 
cording to this model SMBHs form at the center of most galaxies, growing, 
throu gh active star formation, with the mass of their bulges (see Cattanea 



2001 



and references therein). Therefore, one should not be surprised to find 
so ample evidence of AGNs with star formation. The TOs are possibly one 
good example. As the mass of the bulge and of the SMBH grow the galaxy 
transforms into an earlier morphological type. Eventually star formation fades 
away, leaving only a narrow line AGN (a LINER), where a SMBH still accretes 
matter but at a lower rate than before. According to our findings, it is very 
possible that all the narrow-line AGNs passed by a phase of higher activity in 
the past, which could suggest that most LINERs were Sy2s sometime in the 
past, and most probably, the Sy2s will eventually end their life as LINERs. 

The standard AGN model explains very well the i^-shapcd distribution 
traced by the NELGs in the BPT-VO diagram. This distribution reflects the 
accumulation of mass at the center of galaxies through the formation of a 
SMBH and a massive bulge. On the left branch of the ^-shaped distribution, 
we find galaxies with small binding energies. These galaxies form most of their 
stars in a disk and do not develop massive bulges. The relatively small-mass 
SMBH that form in their nucleus do not accrete much matter either, which 
explains why we do not detect them. Our own galaxy is possibly of this kind. 
On the right branch of the z/-shaped distribution, going from left to right, we 
find galaxies with increasing binding energies, and consequently more massive 
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bulges and actively accreting SMBHs. The TOs, Sy2s and LINERs in this 
scenario may be related to different mass scales in the formation of a SMBH. 
Alternatively they may represent different evolutionary phases which would 
explain why the exact boundaries between the activity types are difficult to 
delineate. The LINERs in particular may represent the end products of this 
evolutionary process. 

Within the standard AGN model interpretation the Sy2s and LINERs 
seem to show clear evidence of evolution (within each type): as the SMBH 
grows in mass its accretion rate decreases. This would suggest that the acti- 
vity phase of the SMBH in narrow-line AGNs is of longer duration than usu- 
ally assumed based on quasa r models, varying o ver a period of time of the 
order of a Gyr or more (e.g., Kewlev et al. 20061) . larger than the time nec- 
essary to form their bulges. In fact, in the narrow-line AGNs we may still 
see SMBHs growing in mass by accretion, in parallel with the host galaxies 
building their bulges through star formation (the Sy2s, and possibly also the 
TOs, are good examples). This longer activity timcscale however may not 
apply to the more luminous AGNs like the quasars, suggesting possibly diffe - 
rcnt accretion regimes (e.g., ISmall fc Blandfordlll992t Hopkins et aDl2006bh . 
In particular, we see no evidence of negative feedback from the act ivity of the 



SMB H, as suggested to regulate the activity of the quasars (see iPage et al 



2012L and reference therein). This model would not explain the variation of 



bulge characteristics observed in the narrow-line AGNs, passing from the TOs 
to the Sy2s and LINERs, which is more consistent with a positive feedback. 
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A. RESULTS OF STATISTICAL TEST 
The statistical framework of the tests we u s ed is b ased on a new parametric 



ANOVA mo del introduced by 



Hothorn et al.1 (|2008l ) and developed for the R 



softwarjfl by Herberich et al. ( 2010l) . The max-t test docs the simultaneous 



pairwise comparisons of means under control of the family-wise error rate, 
which is the probability of falsely rejecting the initial hypothesis (i.e., finding 
a significant difference among the means of any two groups in the data set even 
though there is actually no difference present). The test takes into account 
possible hctcrosccdasticity and the unequal sizes of the groups. 

We present the results of the tests in the form of simultaneous confidence 
intervals for all pairwise comparisons of group means. Confidence intervals 
including zero indicate no statistically significant differences. Confidence in- 
tervals near zero suggests some level of similarity, and the farther from zero 
they are the more significant the differences. Note the smallness of some of 
the confidence intervals. This is due to the very large number of data used in 
each bin, which makes the tests extremely significant. 

In Figure[17]we show the confidence intervals associated with Figured] (Sec- 
tion[5]) where we compare the variations of physical parameters in NELGs with 
different activity types. In a) we see the VMA increasing in the sense, SFGls 
=> SFG2s => TOs Sy2s => LINERs, with a significant minor difference be- 
tween the TOs and Sy2s; b) the FWHM increases in the same direction, with 
a minor difference between the Sy2s and LINERs. 

In Figure [15] we show the confidence intervals associated with Figure [T^] 
(Section [5]), where we compare in a) the masses of the SMBHs and in b) the 
accretion rates. We see that LINERs have significanlty more massive SMBHs 
that accrete at lower rates than the Sy2. 



5 http://CRAN.R-project.org 
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Fig. 17. The confidence intervals associated with Figure [5] in Section [5] a) VMA; b) 
FWHM of the Ha. 
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Fig. 18. The confidence intervals associated with Figure [12] in Section [6] a) the 
masses of the SMBHs; b) the accretion rates. 



